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We have studied a simple effective model of charge ordered insulators. The tight binding Hamil- 
tonian consists of the effective on-site interaction U and the intersite density-density interaction 
Wij (both: nearest-neighbor and next-nearest-neighbor) . In the analysis of the phase diagrams and 
thermodynamic properties of this model we have adopted the variational approach, which treats the 
on-site interaction term exactly and the intersite interactions within the mean-field approximation. 
Our investigations of the general case (as a function of the electron concentration n) have shown 
that the system exhibits various critical behaviors including among others bicritical, tricritical, 
critical-end and isolated critical points. In this report we concentrate on the metastable phases and 
transitions between them. One finds that the first- and second order transitions between metastable 
phases can exist in the system. These transitions occur in the neighborhood of first as well as second 
order transitions between stable phases. For the case of on-site attraction the regions of metastable 
homogeneous phases occurrence inside the ranges of phase separated states stability have been also 
determined. 

PACS numbers: 71.10.Fd, 71.45.Lr, 64.60.My, 64.75.Gh, 71.10.Hf 



I. INTRODUCTION 

There is intense research in the field of electron 
charge ordcrings phenomena due to their relevance for 
a broad range of important materials such as mangan- 
ites, cuprates, magnetite, several nickel, vanadium and 
cobalt oxides, heavy fermion systems and numerous or- 
ganic compounds (Refs. [1-5] and references therein). 

The effective Hamiltonian of an electron system on the 
lattice in the zero bandwidth limit considered in this re- 
port can be written in the following form: 
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where cf a denotes the creation operator of an electron 
with spin a at the site i, hi = ^2 a h ia , h ia — cf a Ci a , 
U is the on-site density interaction, W\ and Wi are 
the intersite density-density interactions between near- 
est neighbors (nn) and next-nearest neighbors (nnn), re- 
spectively These interactions will be treated as the effec- 
tive ones and will be assumed to include all the possible 
contributions and renormalizations. fi is the chemical 
potential, depending on the concentration of electrons 
n = J2i ("») ! w i tn < n < 2 and N is the total num- 
ber of lattice sites. Our denotations: tiq = \(ua —tib), 
n a = jj Sj£a(^i); an d a ~ A, B labels the sublattices. 
W = z\W\ + Z2W2, Wq — —Z1W1 + z 2 W 2 , where z\ 
and Z2 are the number of nn and nnn, respectively. 
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We have performed extensive study of the phase dia- 
grams of the model (1) for W\ > and arbitrary n [6-10]. 
Depending on the values of model parameters the system 
can exhibit not only several homogeneous charge ordered 
(CO) phases and nonordered (NO) phase, but also vari- 
ous phase separated (PS) states (PS1: CO-NO, PS2: CO- 
CO, PS3: NO-NO) [8-11], in which two domains with 
different concentration exist (coexistence of two homo- 
geneous phases). However, the behaviors of metastable 
phases occurring in model (1) have not been analyzed till 
now. 

In the analysis we have adopted a variational ap- 
proach (VA), which treats the on-site interaction term 
(U) exactly and the intersite interactions (Wij) within 
the mean- field approximation (MFA). One obtains two 
equations for n and jiq , which are solved self-consistcntly. 
Explicit forms of equations for the free energy and other 
thermodynamical properties are derived in Ref. [8]. uq 
is non-zero in the charge-ordered phase, whereas in the 
nonordered phase nq = 0. Only the two-sublattice or- 
derings on the alternate lattices are considered in this 
report. 

In present report we will concentrate on the possibility 
of metastable phases occurrence on the phase diagrams 
of model considered. 



II. RESULTS AND DISCUSSION 

A. Wi > 0, W 2 = 

For W\ > and W2 = the system exhibits a tricrit- 
ical line, a critical end point line and a line of isolated 
critical points [7]. The CO-NO transition can be second 
order as well as first order. Two different CO phases (i.e. 
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FIG. 1. Phase diagrams for Wi > 0, W 2 = and: (a) U/W = 0.8, (b) U/Wo = 1.0, (c) n = 1. Dotted and solid lines denote first 
and second order transitions between stable phases. Dashed-dotted lines denote the boundaries of metastable phase occurrence 
(names of metastable phases in brackets). Dashed line (panel (a)) denotes second order boundary between metastable phases. 




FIG. 2. Phase diagrams for U/(-Wq) = -10, Wi > and k = z 2 W 2 /zi Wi = -0.2, -0.6, -1.0 (as labeled). CO denotes now 
the HCO phase. Solid and dashed lines indicate second order and "third order" transitions, respectively. Below dashed-dotted 
lines all homogeneous phases are unstable. Details in text. 



LCO and HCO) are separated by first order line. 

In Fig. 1 we presents a few particular phase diagrams 
involving metastable phases. It is quite obvious that 
metastable phases are present in the neighborhood of 
first order (HCO-LCO and HCO-NO) transitions (such 
region is very narrow for the HCO-LCO transition). 
Above the first order transition temperature the phase, 
which was stable below the transition temperature, is 
metastable, and inversely, below the transition temper- 
ature the phase, which was stable above the transition 
temperature, is metastable. However, one should no- 
tice that second order LCO-NO transition occurs be- 
tween two metastable phases with increasing tempera- 
ture connected with continuous change of charge-order 
parameter in metastable phases (Fig. la, U/Wo = 0.8). 
Such transition between metastable phases occurs in the 
higher energy branch of solutions, whereas the lowest en- 
ergy solution is the HCO phase. Other interesting fea- 
ture of the model is that in the vicinity of second or- 
der LCO-NO transition for n > 0.8 the HCO phase is 
metastable (Fig. lb, U/Wq = 1). This behavior is con- 
nected with HCO-NO transition occurring for U/Wo < 1 
(cf. Fig. lc). 

Let us stress that we found all MFA solutions of the 
model considered. Thus metastable phases occur only in 
the regions explicitly denoted on the phase diagrams. In 
other regions there are no metastable phases - only one 



(stable) solution exists. 

B. Wi > 0, W 2 < 

In such range of model parameters the system can ex- 
hibits not only several CO phases, but also various phase 
separated states: PS1 and PS2 [8, 9]. Examples of the 
ksT vs. n phase diagrams evaluated for strong on-site at- 
traction U/(—Wq) = —10, W\ > and various ratios of 
k = Z2W2I Z\W\ < are shown in Fig. 2. A transition 
between homogeneous phase and PS state is symboli- 
cally named as a "third order" transition. During this 
transition a size of one domain in the PS state decreases 
continuously to zero at the transition temperature. The 
CO and NO phases are separated by the second order 
transition line and for k = no metastable phases occur. 
If k < in the ranges of PS stability the homogeneous 
phases can be metastable (if dfi/dn > 0) as well as un- 
stable (if dnjdn < 0). 

For k < —0.6 the PS1 state occurs on the phase dia- 
gram and the critical point for the phase separation (de- 
noted as T) lies on the second order line CO-NO. As 
k — > —00 the T-point occurs at n — 1 and the homoge- 
neous CO phase does not exist beyond half-filling. If 
k = —0.6 iJ-point is present on the phase diagram and 
the system changes a tricritical behavior (for k < —0.6) 
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FIG. 3. Phase diagram for U/(-W Q ) = -10, Wi < and 
W2 = 0. Dashed line indicates the PS3-NO transitions. Be- 
low dashed-dotted line the homogeneous NO phase is unsta- 
ble. 



into a bicritical behavior (for > k > —0.6). In the 
ranges of PS1 stability the NO phase (in region X) 
and the CO phase (in region Y) are metastable. Below 
dashed-dotted lines all homogeneous phases considered 
(CO as well as NO) are unstable (i.e. dfi/dn < in all 
homogeneous solutions) . 

When —0.6 < k < a transition between PS state and 
homogeneous phase takes place at low temperatures, 
leading first to phase separation into two coexisting CO 
phases (PS2), while at still lower temperatures CO and 
NO phases coexist (PS1). The critical point (denoted 
as B) for this phase separation is located inside the CO 
phase. The E-F solid line is associated with continuous 
transition between two different PS states (PS1-PS2, the 
second order CO-NO transition occurs in the domain 
with lower concentration). Similarly as for k < —0.6, 
for —0.6 < k < the NO phase (in region X) or the 
CO phase (in regions Y and Y') are metastable in the 
ranges of PS1 stability. One should notice that second 
order transition CO-NO between metastable phases oc- 
curs (the solid line between regions X and Y' in Fig. 2 
for k — —0.2). At higher temperatures, in the ranges of 
PS2 stability only the CO phase can be metastable (in 
regions Z and Z'). Below dashed-dotted line all homo- 
geneous phases considered are unstable. 

For larger values of U/{— Wq) (especially if 
U/(Wq) > it could be possible that more than 
one metastable phase exist in ranges of PS states 



occurrence, however we do not analyze it in this report. 



C. Wi < 0, W 2 = 

For the case Wi < (W 2 = 0) the model (1) (W 2 = 0) 
exhibits a phase separation NO-NO (electron droplets 
state - PS3) at low temperatures [11]. In this PS state 
different spatial non-ordered regions have different aver- 
age electron concentrations. In such a case, at higher 
temperatures only the homogeneous NO phase occurs. 
The phase diagram for £//|Ho| = — 10 and W\ < in- 
volving metastable phases is shown in Fig. 3. One can 
notice that the homogeneous NO phase is metastable in 
regions V and V . The line restricting (meta-)stability 
of the NO is tangent to the PS3-NO boundary in the R- 
point (R is a bicritical point). Below dashed-dotted line 
the homogeneous NO phase is unstable (i.e. d/i/dn < 
in the NO phase). 



III. CONCLUSIONS 

In this report, we have presented some particular phase 
diagrams of the extended Hubbard model with intersite 
density-density interactions in the zero-bandwidth limit. 
We have found that the first- and second order transi- 
tions between metastable phases can exist in the system. 
These transitions occur in the neighborhood of first as 
well as second order transition between stable phases. 
We have also determined the regions of metastable ho- 
mogeneous phases occurrence inside the ranges of phase 
separated states stability for the case of on-site attrac- 
tion. 
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